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A  MULTIVARIATE  TEST  TO  CHECK  THE  VALIDITY 
OF  A  CERTAIN  STRUCTURE  OF  THE  COVARIANCE 
MATRIX  IN  SCHEFFE' 'S  MIXED  MODEL 

A.  I.  KHURI* 

ABSTRACT 

In  this  article,  Scheffer's  mixed  two-way  classification  model  for 
balanced  as  well  as  unbalanced  data  is  considered.  The  covariance 
structure  associated  with  this  model's  random  effects  is  quite 
general.  Under  the  usual  assumptions  of  independence  of  the 
random  effects  and  homogeneity  of  their  respective  variances,  the 
covariance  structure  assumes  a  particular  form.  A  multivariate 
test  is  presented  to  check  the  validity  of  this  form.  A  numerical 
example  is  given  to  illustrate  the  implementation  of  the  proposed 
test. 

KEY  WORDS:  Balanced  data;  Unbalanced  data;  Mixed  two-way 
classification  model;  Model  I;  Model  II;  Random  effects. 
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1 .  INTRODUCTION 

Different  sets  of  assumptions  can  be  found  in  the  statistical 
literature  concerning  the  random  effects  in  a  mixed  linear 
model.  A  discussion  of  these  assumptions  with  regard  to  the 
balanced  two-way  classification  model  was  given  by  Hocking 
(1973).  the  least  restrictive  of  these  sets  of  assumptions  is  the 
one  considered  by  Scheffe-*  (1956;  1959,  Chapter  8).  In  another 
set  of  assumptions  that  is  commonly  found  in  the  literature,  the 
random  effects  are  described  as  being  statistically  independent 
and  normally  distributed  with  zero  means  and  variances  that  are 
equal  for  all  levels  associated  with  a  particular  random  effect. 

We  shall  refer  to  the  mixed  two-way  classification  model  as  Model 
IX  or  as  Model  I  depending  on  whether  the  latter  assumptions  or 
those  considered  by  Scheffe'  are  valid,  respectively. 

Let  us  consider  the  unbalanced  mixed  two-way  classification 

model 


7ijk  "  u  +  °i  +  +  (®S)ij  +  eijk’  (1*0 

i*l ,2 , . • • ,r ;  j“l ,2, • • » ,s ;  k*l ,2, . . . ,n^j 

where  u  and  are  fixed  unknown  parameters;  8j,  (aS)^,  and 
are  random  variables  assumed  to  be  normally  distributed  with  zero 
means.  For  reasons  to  be  seen  later,  it  is  also  assumed  that  n.. 
>  2rs,  where  n..  -  I ^  ^n^  is  the  total  number  of  observations. 


Let  yij 


be  the  cell  sample  mean  (i*l ,  2 , . . .  ,r; 

,s).  From  (1.1)  we  have 


U  +  +  8^  +  (aB)^  +  eijt» 


i“l ,2, . . . ,r; 
j  =  l ,2,  •  • •  ,s  , 


where 


ij 

( 1  /  n  )  l  t 


ij 


k=  1 


ij^* 


(1.2) 


Model  (1.2)  can  be  written  in  vector  form  as 


-Kir+8  +  Sj  +ly  i'1-2 . .  (K3> 

where 

iS  '  (5U.’  J2J. . J’1,Z . . 

1  is  a  column  vector  of  ones  of  order  r* 1 , 

^r 

a  -  a2 . Oj.)'* 

Sj  “  ® j4>r  +  j  » 

where 

(ciB)j  *  [ (aB) ^  ,  (afl^j  »  •••*  (a6)rj]  *  jml»2,...,s. 


e . 


j-1,2. 


s. 


Under  Model  I,  the  gj  and  fj(j-l ,2, . . . ,s)  are  independently 
distributed  as  N(0,  £)  and  n(Q,  ) ,  respectively,  whre  £  is  the 
variance-covariance  matrix  of  ,  o|  is  the  error  variance,  and 
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is  the  diagonal  matrix 


Kj  “  diag(n‘j,  n^1 ,  ...»  n^1), 


,s. 


(1.4) 


Under  Model  II,  which  is  a  special  case  of  Model  I,  E  has  the  form 


-  erf  J  + 
0  ~r 


j2  I  » 
a8~r 


(1.5) 


where  and  o2g  are  the  variances  of  the  Bj  and  the  (a0)^, 
respectively,  and  Jr  and  ,Jr  are,  respectively,  the  matrix  of  ones 
and  the  identity  matrix,  both  of  order  rxr. 

Quite  often,  when  Model  II 's  assumptions  are  stated  in  the 
literature,  no  attempt  is  made  to  check  their  validi-ty  for  the 
data  under  cosideration.  As  was  mentioned  earlier,  the  difference 
between  Model  I  and  Model  II  is  in  the  structure  of  the  variance- 
covariance  matrix,  £.  If  this  structure  is  of  the  form  given  in 
(1.5),  then  Model  II  is  considered  valid  provided,  of  course,  that 
the  other  assumptions  of  Model  I  are  also  valid.  It  is  Important 
here  to  note  that  Model  II  improses  severe  restrictions  on  the 
variances  and  covariances  of  the  random  effects.  These 
restrictions  may  not  be  easily  justified  in  practice  (see  Scheffe' 
1959,  pp.  264-265;  Alalouf  1980,  p.  195).  The  main  attractive 
feature  of  these  restrictions,  particularly  in  the  balanced  case, 
lies  in  their  legitimization  of  the  use  of  the  usual  F-ratios  for 
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the  traditional  analysis  of  data*  Furthermore,  these  restrictions 
are  superfluous  since  in  the  balanced  case  at  least,  other  exact 
tests  exist  for  the  more  flexible  Model  I  as  was  demonstrated  by 
Scheffe"  (1959,  Chapter  8).  The  latter  tests  are  not,  in  general, 
based  on  using  the  F-ratios  in  the  corresponding  analysis  of 
variance  (ANOVA)  table.  For  example,  the  test  concerning  the 
fixed  effects  was  developed  on  the  basis  of  Hotelling  T2- 
s tat is  tic.  The  tests  concerning  the  random  effects,  however,  were 
obtained  in  the  usual  manner  using  the  ANOVA  F-ratios.  But, 
unlike  Model  II,  the  power  of  the  test  for  the  interaction  is  not 
expressible  in  terms  of  the  central  or  noncentral  F-distribution. 

The  purpose  of  this  article  is  to  present  a  test  to  check  the 
validity  of  the  covariance  structure  in  (1.5)  for  the  unbalanced 
model  in  (1.1).  More  specifically,  under  the  assumptions  of 
Scheffe“'s  Model  I,  the  following  hypothesis  will  be  tested: 

H_  :  Z  ™  aJ  +  bl  ,  (1.6) 

0  ~  ~r  ~r 

where  a  and  b  are  nonnegative  constants.  In  particular,  in  case 
of  equal  subclass  frequencies  (that  is,  when  the  n^j  are  equal), 
the  test  provides  a  check  on  Model  II  for  a  balanced  data  set. 

The  proposed  test,  which  is  multivariate  in  nature,  represents  a 
further  demonstration  of  the  utility  of  the  multivariate  approach 
to  mixed  linear  models.  This  approach  was  discussed  in  the 
balanced  case  by  Alalouf  (1980). 


.  TESTING  THE  VALIDITY  OE  THE  COVARIANCE  STRUCTURE  IN  (1.6) 
Let  us  again  consider  the  model  given  in  (1.3).  Model  I  will 
be  assumed,  that  is,  the  0j  and  are  independently  distributed 
as  N(2,  £)  and  N(Q,  respectively,  where  is  the  diagonal 

matrix  defined  in  (1.4).  It  is  also  assumed  that 

n..  >  2rs,  (2.1) 

where  n..  is  the  toal  number  of  observations,  r  and  s  are  the 
numbers  of  levels  associated  with  the  main  effects. 


The  mean  vector  and  variance-covariance  matrix  of 


in  (1.3) 


+  2>  3ml>2 . .  <2*2> 

Var(^  )  -  Z  +  ,  J-1,2 . .  (2.3) 

The  error  sum  of  squares  for  model  (l.l)  is  given  by 


Q 


(2.4) 


..Formula  (2.4)  can  be  rewritten  as 

Q  -  X 5Z>  (2*5) 

where  ^  ts  the  vector  of  a11  observations 

with  being  the  vector  of  observations  obtained  under 
j (“1 »2, * . * ,s) •  The  matrix  £  is  (2.5)  is  idempotent  of  order 
n..xn«.  and  rank  n..-rs,  and  has  the  form 
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where  J  is  Che  matrix  of  ones  of  order  n. .xn. .  (i=l , 2, . . .  ,r ; 

~n  lj  ij 

ij 

j=l ,2, . . . ,s) .  The  second  term  on  the  right-hand  side  of  (2.6)  is 

the  direct  sum  of  the  J  /n...  Since  R  is  symmetric  with 

~n. .  ij  ~ 

ij 

eigenvalues  equal  to  one  or  zero,  it  can  be  expressed  as 


R  =  C  diag( I  ,  0)  C' , 

~  ~  v~n..-rs  ~ 


(2.7) 


where  "diag"  denotes  a  block-diagonal  matrix,  0  is  a  zero  matrix 
of  order  rsxrs,  and  ,£  is  an  orthogonal  matrix  of  order  n..xn..  . 
If  condition  (2.1)  is  satisfied,  then  (2.7)  can  be  rewritten  as 


R  *  C  diag( I  ,  I 

~  v,'vrc  ^ 


rs  ~n#*-2rs 


•  °)£' 


(2.8) 


The  matrix  £  in  (2.8)  can  be  partitioned  as 


C=[C  :  C  :  ...  :C  tC  tC  1  , 
~  ~1  ~2  ~s  ~s+l  ~s+2 J 


(2.9) 


where  is  a  matrix  of  order  n..xr  for  i»l,2,...,s;  n. .x (n. .-2rs) 
for  i-s+1;  and  n..xrs  for  i»s+2. 

Let  us  now  define  the  vector  w.  as 


■  i,  +  0ir  -  . . 


(2.10) 
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where  y  and  K,  are  given  in  (1.3)  and  (1.4),  respectively,  X  is 

*j  ~j 

given  by 

X  »  l/(min(n  )),  (2.11) 

i.j  J 

and  (XIf  -  K.)*5  is  a  diagonal  matrix  whose  ith  diagonal  element  is 
(X  -  »  i“l»2,...,r;  j*l,2,...,s. 

The  following  theorem  gives  distributional  properties  of  the 
w.  in  (2.10): 

Theorem  2.1.  Under  Model  I,  the  w^  defined  in  (2.10)  have  the 
following  properties: 

(i)  They  are  statistically  indpendent  and  distributed  as 
normal  random  vectors. 

(ii)  E(w  )  -  n  ,  j"l,2, ... ,s. 

(iii)  Var(w  )  -T,  j*l,2,...,s, 

where 


n  -  pi  +  a,  (2.12) 

r  -  Z  +  Xa2 I  .  (2.13) 

<v  e~r 

Proof.  The  proof  is  given  in  the  Appendix. 

It  is  easy  to  verify  that  £  has  the  form  described  in  (1.6)  if  and 
only  if  £  »  aJ  +  (b+Xo2)l  ,  which  is  of  the  form 
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r  =  a'J  +  b'l  , 
~  ~r  ~r 


(2.14) 


where  a'  and  b'  are  nonnegative  constant  coefficients.  Formula 
(2.14)  can  be  rewritten  in  the  more  familiar  form 


r 


=  a 2 


1  P  P 
P  1 
P 


P 

P 


(2.15) 


P 


P 


P 

1 


where  o^-  =  a'  +  b'  and  p  =  a' /(a'  +  b").  We  note  that  (2.15)  is  a 
special  case  of  the  more  general  expression 


r 

«■%/ 


(2.16) 


where  the  correlation  coefficient,  is  not  necessarily 

nonnegative  as  p  is  in  (2.15). 

According  to  Theorem  2.1,  the  rx 1  vectors,  w  ,  w„ ,  ...,  w 

~  1  ~  2  ~  s 

form  a  random  smaple  of  size  s  from  N(n,  £).  This  sample  will  be 
utilized  to  test  the  hypothesis 


The  test  is  based  on  the  likelihood  ratio  principle  and  the  test 
statistic,  given  by  Wilks  (1946,  page  260)  and  denoted  by  L,  is 
such  that 


(2.18) 


(see  also  Seber  1984,  p.  95).  The  rxr  matrix,  £,  in  (2.18)  is  an 
unbiased  estimator  of  £  given  by 

£  ■  2  tej  -  SHxj  -  5)',  (2-i9) 

j«i  j  j 

where 

1  3 

2  *  7  Is-  (2.20) 

j-1 J 

In  (2.18),  |*|  and  tr(»)  denote  the  determinant  and  the  trace  of  a 
matrix.  The  null  hypothesis  in  (2.17)  is  rejected  at  the  a-level 
of  significance  if 

M  <  Ta(r,s), 


(2.21) 


where  x  (r,s)  is  the  lower  100a%  point  of  the  null  distribution  of 
the  statistic  M  in  (2.18).  The  exact  null  distribution  of  M  was 
determined  by  Wilks  (1946)  for  r*2  and  r=3.  He  showed  that  the 
probability  integral  of  M  for  r=2  and  that  of  /M  for  r=*3  are 
incomplete  beta  functions  with  parameters  p=(s-2)/2,  q=l/2  and 
p=s-3,  q=2,  respectively.  On  that  basis,  he  obtained  the 
percentage  points  of  M  for  r*2  and  r=3  using  the  beta  distribution 
tables  (see  Wilks  1946,  Table  1,  p.  263).  Nagarsenker  (1975) 
derived  the  exact  distribution  of  M  and  provided  a  tabulation  of 
the  corresponding  5%  and  1%  points  for  r»4(l)10  and  various  values 
of  s.  For  large  values  of  s  and  under  the  null  hypothesis  (2.17), 
the  statistic, 

V  =  -{ s-l-[ r( r+1 )2 (2r-3) ] / [ 6 ( r— 1 ) ( r2+r-4)] } log(M) ,  (2.22) 

is  asymptotically  distributed  as  chi-square  with  [r(r+l)-4]/2 
degrees  of  freedom  (see  Box  1950,  p.  375;  Seber  1984,  p.  95).  In 
(2.22),  log(M)  denotes  the  natural  logarithm  of  M. 

In  particular,  if  the  data  set  is  balanced,  that  is,  when  njj 
=*  n  for  all  i,j,  then  from  (1.4)  and  (2.11)  we  have 

\Ir  “  “  0»  j“l»2,...,s. 

Hence,  w,  -  y.  for  j-l,2,...,s  in  (2.10).  In  this  case,  the  test 
~J  J 

statistic  in  (2.18)  can  be  used  with  ^  substituted  for 

Jgj(j“l*2 . . 


Rejection  of  the  null  hypothesis  in  (2.17)  gives  support  to 

the  contention  that  £  is  not  of  the  form  described  in  (2.16), 

hence,  it  cannot  have  the  form  given  in  (2.15).  This  in  turn 
implies  that  £  is  not  expressible  as  in  (1.6).  Consequently, 

rejection  of  the  null  hypothesis,  Hq,  in  (2.17)  is  an  indication 

that  the  null  hypothesis,  Hq,  in  (1.6)  is  false,  which  can  be 
construed  as  evidence  that  Model  II  is  invalid.  Failing  to  reject 
the  null  hypothesis  in  (2.17),  however,  does  not  necessarily  mean 
that  Model  II  is  valid.  This  is  because  this  hypothesis  does  not 
imply  the  one  in  (1.6).  In  the  latter  hypothesis,  both  a  and  b 
must  be  nonnegative.  From  (2.17)  and  (2.13)  we  can  only  conclude 
(1.6),  but  without  the  nonnegativity  constraint  on  both  a  and  b. 

3.  A  NUMERICAL  EXAMPLE 

We  consider  an  example  given  by  Raid  (1952,  p.  471)  and 
reproduced  in  Scheffe'  (1959,  p.  289,  problem  8.2).  In  this 
example,  measurements  were  given  of  the  waterproof  quality  of 
sheets  of  material  manufactured  by  r-3  different  machines  (factor 
A)  over  a  period  of  s-9  days  (factor  B).  Three  sheets  were 
selected  for  each  machine-day  combination.  The  response  values, 
namely,  the  logarithms  of  the  permeability  in  seconds  for  the 
sheets,  are  given  in  Table  l.  We  refer  to  this  balanced  data  set 
as  Data  Set  1.  The  effects  associated  with  factors  A  and  B  are 
considered  fixed  and  random,  respectively. 
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An  unbalanced  data  set  (Data  Set  2)  was  obtained  by  deleting 
certain  response  values  from  Data  Set  1.  The  deleted  values  are 
marked  with  an  asterisk  in  Table  1.  The  matrix  r  in  (2.19)  for 
Data  Set  2  has  the  value 


r  = 


.0239  .0038  .0072 
.0038  .0050  -.0038 
.0072  -.0038  .0160 


From  (2.18)  we  find  that  M  =  L2^9  =  .2806.  From  Beyer's  (1968) 
tables  of  percentage  points  of  the  beta  distribution  (Table 
III. 10,  p.  256),  the  corresponding  lower  10%  point  of  the 
statistic  M  is  equal  to  (.54744)2  =  .2997  (recall  that  for  r*3,  /M 
has  the  beta  distribution  with  parameters  p  *  s-3  =  6  and  q  = 

2).  The  test  is,  therefore,  significant  at  the  10%  level.  We 
thus  conclude  that  Model  II  is  invalid  for  Data  Set  2. 

Let  us  consider  next  Data  Set  1  (this  is  the  original 
balanced  data  set  in  Table  1).  In  this  case. 


r 


.0288 

.0034 

.0012 


.0034 

.0067 

.0005 


.0012 

.0005 

.0236 


and  M  ■  .5737,  which  exceeds  the  10%  critical  value,  .2997. 


Consequently,  we  cannot  conclude  invalidity  of  Model  II  for  Data 


Set  1  (this  does  not  mean  that  Model  II  is  valid). 

We  now  show  that  a  slight  modification  in  Data  Set  1  can  lead 
to  a  highly  significant  test,  we  shall  keep  the  same  subclass 
frequencies  as  in  Data  Set  1,  but.  alter  only  seven  response  values 
in  Table  1.  The  new  response  values  are  given  in  parentheses  next 
to  the  values  that  have  been  altered.  We  note  that  the  new  values 
do  not  differ  by  much  from  their  original  counterparts.  The  new 
balanced  data  set  thus  obtained  is  called  Data  Set  3.  The 
matrix  r  now  has  the  value 


.0300  .0105  -.0049 

£  -  .0105  .0083  -.0032 

-.0049  -.0032  .0025 

The  corresponding  value  of  M  in  (2.18)  is  .0677.  This  is  a  highly 
significant  result  since  the  lower  .5Z  point  of  the  statistic  M 
from  Beyer's  (1968,  p.  264)  tables  is  (.31509)2  .  .0993.  Thus,  a 

minor  change  in  Data  Set  1  has  led  to  a  rejection  of  Model  II.  We 

can,  therefore,  clearly  see  the  Importance  of  checking  the 
validity  of  Model  II's  assumptions  concerning  £  prior  to  using  the 

standard  techniques  of  analysis  of  variance  for  data  analysis. 


APPENDIX:  PROOF  OF  THEOREM  2.1 


(i)  It  is  obvious  that  the  w^  in  (2.10)  are  normally 
distributed  ( j-1 , 2, . . . ,s) .  Now,  the  two-way  classification  model 
in  (l.l)  can  be  written  as 

y  -  Ul  +  X.o  +  X.8  +  X,(a$)  +  e,  (A. 1 ) 

~  ~n .  •  ~1~  ~2~  ~3  ~  ~ 

where  X  ,  X  ,  and  X_  are  matrices  of  zeros  and  ones  of  orders 

I  t.  ^ 

n..  x  r,  n..  x  s,  and  n..  x  rs,  respectively;  g,  £ ,  and  (a8  )  are 
vectors  whose  elements  are  the  c^,  8 j ,  and  (aS)^,  respectively. 
The  variance-covariance  matrix  of  ^  is 

i  •  [*2  :  *31  Var<*)(*2  •  It,]**  (*-2) 

where  J  ■  [s',  (as)']'.  Let  %  *  (Xl»  22»  ***’  1)6  t*le  veccor 

of  cell  sample  means,  where  £  (j-l,2 . s)  is  defined  in 

(1.3).  This  vector  can  be  expressed  as 

(A. 3) 


where  g  is  the  direct  sum 
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Let  us  now  consider  Che  matrix  &  defined  in  (2.6).  From  (2.8)  and 
(2.9)  we  can  write 


s+1 

*  -  I  Si  s:  .  (a.5) 

j-i  3  3 


where 


C'  C,  -  I  ,  j*l ,2, . . . ,s+l 
~j  ~j  ~r 

C'  -  0,  j*j'. 

-yj  >v  *  w  J 


(A.6) 


It  is  easy  to  verify  that 


DR  -  0,  RX.  -  0,  RX,  -  0. 


(A.  7) 


From  (A. 2)  and  (A. 7)  we  conclude  that 


D  A  R  -  0. 

y  m 


(A. 8) 


Also,  from  (A.5),  (A.6),  and  (A. 8)  we  have 


D  A  C  -  0, 
—  — 


, 2 , . . . ,s. 


(A. 9) 


Moreover, 
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Cov(i»  yCj)  *  D  A  Cj  ,  j-l,2,...,s. 


(A. 10) 


Formulas  (A. 9)  and  (A. 10)  imply  that  %  and  Cl  and 

hence  v  and  C'  are  statistically  independent  for  2., 

j«l,2,...,s.  Furthermore,  under  Model  I,  the  jr  ( j  =  l , 2 , . . .  ,s)  are 
statistically  independent.  The  C'  %  (j=*l  ,2, . . .  ,s)  are  also 
statistically  independent  due  to  the  fact  that 


Cov(£j  X  Sr)  *  A  Cy  =  0,  j*r 


( A. 1 1 ) 


since 


CJ  *2  "  2-  Sj  h  *  2*  (A. 12) 

SJ£r-2»  > 

by  (A. 5),  (A. 6),  and  (A. 7).  From  the  above  arguments  we  conclude 
that  the  ( j-1 , 2 , . . . , s)  in  (2.10)  are  statistically  independent, 
(ii)  From  (A.l)  we  have 


E(Wj)  -  E(X)M\Ir  -  +  Xja). 


(A. 13) 


Now,  in  addition  to  (A. 7),  the  matrix  £  satisfies 


R  1 

~  ~n. . 


2. 


RX 


1 


0. 


(A. 14) 
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From  (A. 5),  (A. 6),  and  (A. 14)  we  thus  have 


£5  in..  “2.  SJi-S.  ^>2 . 


(A. 15) 


From  (2.2),  (A. 13),  and  (A. 15)  we  conclude  that 

*(Zj)  -<*,)- »ir  ♦  2.  . . 

(ill)  In  part  (i)  it  was  shown  that  ^  and  j  are 
statistically  independent  for  j-l,2,...,s.  It  follows  that 


v«(sJ)  -  +'(»!,  -  «  Sj(”r  -  £j)\ 


j"l ,2,..* ,3. 


(A. 16) 


But,  from  (A.6)  and  (A. 12)  we  note  that 


C 1  A  C. 


Hence,  for  j»l,2,...,s. 


cj2!  , 
e~r* 


jsl,2,...,s< 


var(2j)  "  Var(Ij)  +  a|(Air  -  &j) 

-  Z  ♦  ojKj  +  a|(X!r  “  Kj ) ,  by  (2.3), 


■  Z  *  \a2I 
—  c~r 


r. 


(A. 17) 
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TABLE  1. 

THE  LOGARITHMS 

OF  THE  PERMEABILITIES 

IN  SECONDS 

(Hald 

1952,  p.  472) 

Day 

Machine 

1 

2 

3 

1.404* 

1.306 

1.932* 

1 

1.346 

1.628 

1.674 

1.618 

1.410* 

1.399(1.10) 

1.447 

1.241* 

1.426(1.126) 

2 

1.569 

1.185 

1.768 

1.820 

1.516 

1.859* 

1.914 

1.506 

1.382 

3 

1.477 

1.575* 

1.690 

1.894 

1.649 

1.361 

1.887 

1.673 

1.721 

4 

1.485 

1.372* 

1.528 

1.392 

1.114 

1.371 

1.772* 

1.227 

1.320 

5 

1.728 

1.397* 

1.489 

1.545(1.245) 

1.531 

1.336(1.691) 

1.665 

1.404 

1.633 

6 

1.539 

1.452* 

1.612 

1.680 

1.627 

1.359 

1.918 

“  1.229*(1.729) 

1.328 

7 

1.931 

1.508 

1.802 

2.129 

1.436 

1.385 

1.845 

1.583 

1.689 

8 

1.790 

1.627 

2.248(1.646) 

1.282 

1.795(1.195) 

1.540 

1.636 

1.703 

9 

1.428 

1.067 

1.370 

1.704 

1.384 

1.839 

Note:  Values  marked  with  an  asterisk  are  deleted  for  Data  Set  2. 
Values  in  parentheses  replace  the  altered  values  to  the  left  for 
Data  Set  3. 
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